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Abstract 

Photo-acoustic tomography is a newly developed hybrid imaging modality 
that combines a high-resolution modality with a high-contrast modality. We 
analyze the reconstruction of diffusion and absorption parameters in an elliptic 
equation and improve an earlier result of Bal and Uhlmann [5] to the partial 
date case. We show that the reconstruction can be uniquely determined by the 
knowledge of 4 internal data based on well-chosen partial boundary conditions. 
Stability of this reconstruction is ensured if a convexity condition is satisfied. 
Similar stability result is obtained without this geometric constraint if 4n 
well-chosen partial boundary conditions are available, where n is the spatial 
dimension. The set of well-chosen boundary measurements is characterized 
by some complex geometric optics (CGO) solutions vanishing on a part of the 
boundary. 

1 Introduction 

Typical medical imaging modalities like computerized tomography (CT), magnetic 
resonance imaging(MRI) and ultrasound imaging(UI) have high resolution but low 
contrast. On the other hand, some modalities, based on optical, elastic, or electrical 
properties of tissues, exhibit a sufficiently high contrast between different types of 
tissues, but involve a highly smoothing measurement operator and are thus typically 
low-resolution. Such examples are optical tomography(OP), electrical impedance 
tomography (EIT) and elastographic imaging(EI). 

Hybrid methods combine high resolution with high contrast, based on physi- 
cal coupling mechanisms. Photo-acoustic tomography(PAT) and thermo-acoustic 
tomography(TAT) are recent hybrid methods based on the photo-acoustic effect 
which couples optical and ultrasonic waves. When a body is exposed to short 
pulse radiation, it absorbs energy and expands thermo-elastically. The expansion 
emits acoustic pulses, which travel to the boundary of the domain of interest where 
they are measured. What distinguishes the two modalities is that in PAT, radia- 
tion is high-frequency radiation (near- infra-red with sub- fim wavelength), while in 
thermo-acoustic, radiation is low-frequency radiation (microwave with wavelengths 
comparable to lm). 

The first step in both PAT and TAT is the reconstruction of the absorbed 
radiation, or deposited energy, from time-dependent boundary measurements of 
acoustic signals. Acoustic signals propagate in fairly homogeneous domains as the 



1 



sound speed is assumed to be known. The reconstruction of the amount of deposited 
energy is therefore quite accurate in practical settings. This step has been studied 
extensively in the mathematical literature, see, e.g. [TJ [7J |H1 HI HH EH HH HS1 fT71 

ii mum m w\- 

The second step in PAT or TAT is called quantitative photo- or thermo-acoustics 
and consists of reconstructing the attenuation and diffusion coefficients from knowl- 
edge of the amount of absorbed radiation. This second step is different in PAT and 
TAT as radiation is typically modeled by transport or diffusion equations in the 
former case and Maxwell's equations in the latter case. The problem of interest in 
the paper is quantitative photo-acoustics(QPAT). 

In [5J, G. Bal and G. Uhlmann studied the uniqueness and stability in the recon- 
struction of the attenuation and diffusion coefficients in QPAT. In this paper, we 
extend Bal and Uhlmann's results to partial boundary illumination conditions and 
prove the uniqueness and stability of the reconstruction. Particularly, we show that 
two coefficients in the diffusion equation are uniquely determined by four well-chosen 
illuminations at part of the boundary of the domain and presented an explicit recon- 
struction procedure. The stability of the reconstruction is established from cither 
four internal data under geometric conditions of strict convexity on the domain of 
interest, or from An well-chosen partial boundary conditions, where n is the spatial 
dimension. 

Mathematically, by the standard Liouville change of variables, the diffusion 
equation is replaced by a Schrodinger equation with unknown potential and with 
internal and partial boundary measurements. By adapting the complex geometri- 
cal optics(CGO) solutions proposed in [19], we are able to obtain uniqueness and 
stability results for the inverse Schrodinger problem. The inverse Liouville change 
of variables concludes the uniqueness and stability of the the diffusive regime. 

The rest of the paper is structured as follows. Section 2 presents the PAT 
problem and our main results. The inverse Schrodinger problem and explicit re- 
construction algorithms arc considered in section 3 and 4. The final results on the 
inverse diffusion problem are proved in section 5. 



2 Inverse diffusion problem and main results 

In photo-acoustic diffusion regime, photon propagation is modeled by the second- 
order elliptic equation, 

-V ■ D(x)Vu + o a (x)u = in f2 

u = g on <9f2, 

where f2 is an open, bounded, connected domain in M. n with C 2 boundary DO., D(x) 
is a diffusion coefficient, a a (x) is an attenuation coefficient, and g is the prescribed 
illumination source on <9f2. In this paper, we will consider the partial data problem, 
i.e., g G C k,a (tt) is supported only on a subset of the boundary. Precisely, let 
xq G K"\ch(f2), where ch(f2) denotes the convex hull of SI. We define the front and 
back sides of dQ with respect to xq by 

dQ± = {x G dQ : ±{x - x) ■ n{x) > 0}, (2) 
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where n(x) is the unit exterior norm at x. Let T be an open subset of <9f2, such 
that dfl + C T. We assume that supp(g) C T. 

The energy deposited by the radiation results thermal expansion, which then 
emits acoustic pressure wave p(t,x). The energy deposited is given by 

d(x) = G(x)a a (x)u in f2, (3) 

where G(x) is the Griineisen coefficient. The acoustic wave is modeled by 

' dup-Ap = in (0,T) x W\ 

d t p\t=o = 0, (4) 
p\t=o = f, 

where c(x) is the speed of the acoustic wave, T > is fixed, and / is the unknow 
initial pressure, which is proportional to the deposited energy. The pressure p(t, x) 
is then measured on [0, T] x <9f2, which allows us to reconstruct / by solving a well- 
posed inverse wave problem, and thus to construct the energy deposited;, see e.g. 
[131 HB1 HOI UllllZ]- We assume this step and also assume G(x) in known. The inverse 
problem of interest in this paper is to reconstruct the coefficients (D(x), o~ a (x)) 
from the partial boundary illuminations gj and the internal measurements dj, for 
1 < j < J, J G N*, with Uj solving 

— V • D(x)Vuj + a a (x)uj = in 

Uj = gj on oil. 

The main purpose of this paper is to prove the uniqueness and the stability of 
the coefficient reconstruction. We define the set of coefficients (D(x),a a (x)) G M 
as 

M = {(D(x),a a (x)) : (VD,a a ) eYx C k+1 (Cl), \\VD\\ Y + \K\\ ck+im < M}, (6) 

Where Y = H n / 2+k+2+e (0) C C k+2 {VL) and M > is fixed. 

We define a subset Cl C f2 to be the complement of a neighborhood of <9S1- in 
f2, i.e., il — tt\(J\f(d£l-)), where J\f(dil-) is a neighborhood of <951_ in f2. 

The main results for the inverse diffusion problem with internal data and partial 
boundary data are as follows, where the measurements g and d are real- valued. 

Theorem 2.1. Let O be an open, bounded, connected domain with C 2 boundary 
dd. LetT and be defined as above. Assume that (D(x),a a (x)) and (D(x),a a (x)) 
are in A4 with D\r — D\r- Let d = (dj) and d — (dj), j — f, . . . ,4, be the in- 
ternal data for coefficients (D(x),a a (x)) and (D(x),a a (x)), respectively and with 
boundary conditions g = (gj), j = f,...,4. Then there is a set of illuminations 
g G (C fc ' Q (9ri)) 4 7 supp(g) C T, for some a > 1/2, such that if d = d, then 
(D(x),a a (x)) = (D(x),a a (x)) in tt. 

To consider the stability of the reconstruction, additional geometric information 
about dft is needed. We impose the following hypothesis. 

Hypothesis 2.2. Let SI and 9S1+ be as above. There exists R < oo such that for 
each y G dfl+ C dtt, we have fi C B y (R), where B y (R) is a ball of radius R that is 
tangent to dfl+ at y. 
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The stability result follows. 

Theorem 2.3. Let k > 3. Let satisfy Hypothesis HH with dQ of class C k+1 . 
Let r and fl be defined as above. Assume that (D(x),a a (x)) and (D(x),a a (x)) 
are in M with D\r = D\r- Let d = (dj) and d = (dj), j = 1,...,4, be the 
internal data for coefficients (D(x), o~ a (x)) and (D(x),a a (x)), respectively and with 
boundary conditions g = (gj), j = 1,...,4. Then there is a set of illuminations 
g G (C fc ' Q (90)) 4 , supp(g) C r, for some a > 1/2 and a constant C such that 

\\ D - -Dllc"- 1 ^) + Iko - O'aUck-^o) < C \\ d - d \\(C k (fl)) 4 - ( 7 ) 

Instead of imposing geometric hypothesis, the stability of the reconstruction 
could also follow from An well-chosen measurements. 

Theorem 2.4. Let k > 2. Let £1 be an arbitrary bounded domain with dQ of 
class C k+1 . Let T and f2 be defined as above. Assume that (D(x),o~ a (x)) and 
(D(x),a a (x)) are in M with D\r = -D|r- Let d = (dj) and d = (dj), j = 1, . . . ,4n, 
be the internal data for coefficients (D(x), a a (x)) and (D(x), a a (x)), respectively and 
with boundary condition g = (gj), j = 1, . . . , 4n. Then there is a set of illuminations 
g E (C k ' a (dtt)) en , supp(g) C T, for some a > 1/2 and a constant C such that 

W D - D \\c k (n) + II "" -°a\\c k (n) - c 'll d - d ll(a fc + 1 (n))*»- ( 8 ) 

The proof of these results are mainly based on the corresponding results for the 
inverse Schrodinger equation. In the following section, we will focus on Schrodingcr 
equation and prove similar results. 



3 Inverse Schrodinger equation with internal data 

Let Q be an open, bounded, connected domain in R™ with smooth boundary <9f2. 
Let dtt+ be defined in @ and let L C dil be an open subset such that dft+ C L. 
Consider the diffusion equation with unknow diffusion coefficient D and unknow 
attenuation coefficient o~ a : 

—V • DVu + a a u = 0, in O, u = g, on 90, 

where supp(g) C L. Using the standard Liouvillc change of variables, v = \jDu 
solves 

Av + qv = 0, (9) 

with 



H Vd d 

The internal data in photo-acoustics is given by 



(10) 



d = a a u = -^=v = jjLV, (11) 



while the new boundary condition is given by V Dg on dQ, assuming D is known 
on <9£1. 
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After relabeling, we have a Schrodinger equation of the form 

Aua + qui = in £1 

3 q 1 (12) 
Uj = gj on 9f2, 

where 1 < j < J, J € N* is the number of illuminations, q is an unknown potential, 
and supp(pj) C T. We assume that the homogeneous problem with gj = admits 
the unique solution Uj = so that A = is not in the spectrum of A + q. We 
assume that q on f2 is the restriction to ft of a function q compactly supported on 
K" and such that q e H n l 2+k+£ (W 1 ) with e > for k > 1. Moreover, we assume 
that the extension is chosen so that 

(13) 

for some constant C = C(£l,k,n) independent of q. That such a constant exists 
may be found e.g. in [24], Chapter VI, Theorem 5. 

We assume that gj e C k ' a {dti) with a > \ and dfl is of class C k+1 so that (O 
admits a unique solution u.j e C k+1 (n), see [15], Theo rem 6.19. The internal data 
are of the form 

dj(x) = fi(x)uj(x), in fi, 1 < j < J. (14) 

Here [i e C fe+1 (f2) verifies < ^ < n( x ) — Mo 

The partial data inverse Schrodinger problem with internal data (PISID) consists 
of reconstructing (g,/z) in fl from knowledge of d = (dj)i<j<j E (C k+1 (il)) J and 
illuminations g = (gj)i<j<j G (C k,a (dQ)) J , while supp(<7j) C T. We will mostly be 
concerned with the case J = 4 and J = An with gj and d, real-valued measurements. 

3.1 Construction of complex geometrical optics solution 

The proof mainly depends on the Complex Geometrical Optics(CGO) solutions. 
Let Po = —h 2 A = J2(hD Xk ) 2 be a differential operator, with h a small parameter 
and D Xk = —id Xk . Let ip <G C°°(0;R), with Vip ^ everywhere. Consider the 
conjugated operator 

n 

e v ' h o P o e~^ h = £(fcD XJb + id *^? = A + lB > ( 15 ) 
fc=i 

where A, B are the formally selfadjoint operators, 

A = (hD) 2 -(V<p) 2 , 

B = Vip-hD + hD- V(/9, 

having the principal symbols 

a = £ 2 -(V<^) 2 , b = 2V(p-£. (16) 

We want the conjugated operator e v ' h oP oe^ ip ^ h to be locally solvable in a semi- 
classical sense, which means its principal symbol satisfies Hormander's condition 

{a,b} = 0, when a = 6 = 0, (17) 

where {a, 6} = J2k(^ik a ^x k b — d Xk ad^ k b) is the Poisson bracket. 
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Definition 3.1. A real smooth function ip on an open set Q is said to be a limiting 
Carleman weight if it has non-vanishing gradient on and if the symbols \16\) 
satisfy the condition (|17p . This is equivalent to say that 

(tp'"Vip, Vip) + (tp"£, £) = when ( 2 = (V</?) 2 and • £ = 0, (18) 

where <p" is the Hessian matrix of <f>. 

Let H s (yt) be the Sobolev space. We denote by i/* d (f2) the semi-classical 
Sobolev space of order 1 on with associated norm 

\Hl lel{u) = \\hVuf + \\u\\ 2 (19) 

and by H^ cl (M. n ) the semi-classical Sobolev space of order s on M" with associated 
norm 

h\\H Ll (mn) = \\(hDyu\\ 2 L2{Rn) = J (l + h 2 ey\u(0\ 2 dt. (20) 

Let P = Pq + h 2 q. Carleman estimates give the following solvablity result (TT1 ITS"] : 

Proposition 3.2. Let be a limiting Carleman weight and q £ L°°(Q). Let < 
s < 1. There exists ho > such that for < h < ho and for every w £ 11*^(0), 
there exists u £ H* cl (Q,) such that 



e^Pe = h\\u\\ H s ci (n) < C\\v\\ H s-i 



Now we start to construct CGO solutions by choosing ip £ C°°(f2) to be a 
limiting Carleman weight and ip £ C°°(0) such that 

(Vip) 2 = (V<^) 2 , Vi>-Vip = Q. (21) 

Then ip is a local solution to the Hamilton- Jacobi problem 

a(a:,VV0 = b(x,Vip) = 0. (22) 

Therefore, 

e -j : (- v>+ i^) Poe j-(- v+ i^) a = [((^ + V ^) 2 - Vip 2 )+i(V<p(hD + V<ip) 

+(hD + Vip)V<p)]a 
= (hL-h 2 A)a, (23) 

where L is the transport operator given by 

L = VipD + DVip + i(VcpD + DV<p). (24) 

There exists a non- vanishing smooth function a £ C°°(f2), see [IQl H9] , such that 

La = 0. (25) 

Assume that q £ L°°(Q) and recall P = h 2 (A + q) = P + h 2 q. Then (23]) implies 
that 

Pei^ + ^a = e^ /h h 2 K, 
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with k = 0(1) in L°° and hence in L 2 . Now Proposition 13.21 implies there exists 
r(x,h) S Hl cl (tt) such that 

eV/hpeU-v+^r = ~h 2 K, and , . 

IklU^cn) = H Wr ll + IMI ^ for somc a 

Hence, 

P( e i^ +i ^(a + r)) = 0, 
i.e., we constructed a solution of the Schrddingcr equation of the form 

e i(-<fi+^)( a + r ). (27) 

In our partial data model, we illuminate on part of the boundary, i.e., supp(g) C 
r. But the CGO solution in (|27|) is non-vanishing everywhere. Fortunately, using 
[19], we can modify the CGO solution to match our model. 

Define T_ = d£l\T. Then, T_ is a strict open subset of <9fJ_. 

Proposition 3.3. We can construct a solution of 

Pu = 0, u\ r _ = (28) 

of the form 

u = ei^+^\a + r)+z (29) 

where tp,ip and a are chosen as above and z = e l hb(x; h) with b a symbol of order 
zero in h and 

Im l(x) = —(p(x) + k{x) 
where k{x) ~ dist(x, <9f2_) in a neighborhood of <9f2_ and b has its support in that 
neighborhood. Moreover, ||r||#o =(D(h), r\gn_ =0, ||(Vv? ■ ^) 3 r||aa + = 0{h?). 

Remark: According to the proof in |19j . supp(z) is arbitrarily close to c9£l in 

n. 

Notice that ip is a limiting Carleman weight, so is — (p. We construct the CGO 
solutions of the form 

u 2 = eM-f +v ^{a 2 + r 2 ) + z 2 . 

In particular, we choose 



(30) 



<f(x) = log | x- x \ and ip(x) = d S n-i [ — —,ui ) , (31) 

\x-x \ 



where x e M"\ch(0) and w e 5" _1 . 

Note that Zj, j = 1,2, are supported only in an arbitrarily neighborhood of 
<9S!_. For the rest of this paper, we will mainly consider the uniqueness and the 
stability of the solutions on the subregion defined by 

O = f2\(supp(zi) U supp(z2))- (32) 

Remark that a neighborhood of the point y G dfl+ n <9il_ = {y e dil | — (y — xq) ■ 
n{y) = 0} is excluded from fi. Thus, there exists a fixed constant r\ > 0, such that 
if x G f2 and y G 917+ n <9f2, then — (x — xo) • n(y) > n. In the following proof, we 
mainly focus on f2, where z\ and z 2 vanish. 
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3.2 Construction of vector fields and uniqueness result 

We assume that we can impose the complex-valued illuminations <?j on <9f2 with 
supp(gj) C r and observe the complex- valued internal data dj, j = 1,2. Note that, 
to make up complex- valued gj and dj, j = 1,2, we need to four real observations. 
dj are of the form dj = fxv,j in Jl, where Uj is the solutions of 

Auj + quj = 0, in Uj — gj, on dtt, j = 1, 2. (33) 

Direct calculation gives us that 

U\AU2 — 1*2 Aui = 0. 

We assume that /i € C fe+1 (0) is bounded above and below by positive constants. 
By substituting Uj = dj/fi, wc obtain that 



or equivalcntly, 
where 



2(diVd 2 - d 2 Vdi) ■ - (diAd 2 - d 2 Adi)/x = 0, 

fid ■ + 7rf^ = 0, (34) 

& := x (x)(d 1 Vd 2 -daVdi) 

1 — /3 • V/i (35) 
7d := ~2 X (' T ^ dlArf2 _ rf 2Adi) = . 

Here, x( x ) is an y smooth known complex- valued function with Ix^)! uniformly 
bounded below by a positive constant on Cl. Note that by assumption on /i, we 
have that f3 d g (C*(f2; C))" and 7 e C fc (fi; C). 

The methodology for the reconstruction of (n,q) will be the same as in [5]: we 
first reconstruct fi using the real part or the imaginary part of (|34j) and the boundary 
condition /1 = d/g on 9f2+. Then we may recover Uj = dj/fi, j = 1,2, explicitly 
and therefore q from the Schrodinger equation Q33p . 

The reconstruction of \i is unique as long as the integral curves of (the real 
part or the imaginary part of) (3d join any interior point x £ £1 to a point on 
dfl+ , denoted as This property of j3 is not guaranteed in general, unless 

the boundary conditions gj, j = 1, 2, are properly chosen. CGO solutions are then 
employed to construct a family of boundary conditions gj, j = 1,2, which could 
produce well-performed vector fields fid- 

To make a distinction from the observed data, we denote the CGO solutions by 
u\ and U2, and denote the internal data and the vector field constructed from CGO 
solutions by d and fid, respectively. 

Proposition 3.4. Let fi be the the normalized vector field, defined by 

P=\fa (36) 
There is a open set of g in C k ' a (dil), with supp(g) C r, such that, for any small e, 

<Ch(l + e) on SI, (37) 



\x Q - xr 



c k (o,c) 



where xq £ R ra \ch(J7) and T is a function of class C fc (f2). Therefore, \34\) admits a 
unique solution fi on Q. 
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Proof. Let U\,U2 be CGO solutions in (j3"0")l . By substituting dj = /j,Uj and x( x ) = 
e~x^ in (f3"5j) . we find /3d, restricted to ft, is given by 

Pd = M 2 ( - + ? 'i)( fl 2 + r 2) + (ai + n)(Va a + Vr 2 ) 

- (a 2 +r 2 )(Vai + Vri)). 

We may then define, on ft, 

$ = ~Pd = Vrv^ + A, (38) 

where T = (ai + ri)(a 2 + r 2 ) and # = + ri)(Va 2 + Vr 2 ) - (a 2 + ra)(Voi + 

Vrx)) < Co/i for some constant Co- Also note that, from (|31[) . V0 = \ x °Z x \i ■ 

We will next choose appropriate boundary conditions g^, j = 1, 2, on <9ft, which 
could lead to small \\/3d — in ft. In particular, recall that iij = on c?Sl_ , for 
some e > 0, we choose gj G C k ' a (dfl) and gj = on <9ft-, such that, 

lift - Wjllcfc.-can) < 6 on 50, j = 1,2. (39) 

Let Uj be the solutions of (|33|) with boundary conditions <?.,■ from (|3"9")l . By elliptic 
regularity, we have that 



i j\\C k + 1 (Q,C) 



<Ce on ft, j = l,2, (40) 



for some positive constant C. Notice that dj = fj,Uj and fi G C +1 (ft), we deduce 
that 

\\ d j - dj\\ G "+^(u t c) < Coe on ft, j = 1, 2. (41) 

Thus, restricting to ft, and (|4"Tj) induce ([3T)) . which indicates, when h,e 
are sufheiently small, (3 is close to a non- vanishing vector —a 2 T = u 2 T , x °Z x \ -z ■ 

\Xq X\ 

Approximately, the integral curves of (3 are rays from any x G ft to xq G ]R"\ch(ft). 
intersecting <9ft + at a point x+(x). Therefore, with fi = d\jgi = <i 2 /<7 2 known on 
dft+, 

/3 ■ + jh = 0, 7 = - ld (42) 

provides a unique reconstruction for [i, so does (|34"|) . More precisely, consider the 
flow ^(i) associated to f3, i.e., the solution to 

e x (t) = p(e x (t)), e x (o) = x£&. 

By the Picard-Lindclof theorem, the above equations admit unique solution 6 X while 
j3 is of class C . Since (3 is non- vanishing, 6 X reaches x+{x) G 9ft+ in a finite time, 
denoted as t + (x), i.e., 

9 x {t + {x)) = x + (x). 
Then by the method of characteristics, (J-(x) is uniquely determined by 

H(x) = n (x + (x))e- ^' +<X) T(«.W)«J» ) (43) 
where ^to = <i/g on <9ft + . This finishes the proof. □ 
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Let us define the set of parameters 

V = {(/*>?) e C fe+1 (^) x iJ* +fc+e (rj); is not an eigenvalue of A + q, 

IMIc"=+i(f2) + Ikll^f + fe+e (f2) - p < °°}- 

The above construction of the vector field allows us to obtain the following unique- 
ness result. 

Theorem 3.5. Let be a bounded, open subset o/R™ with boundary of class C k+1 . 
Let (fJ.,q) and (fi,q) be two elements in V and £1 be defined in (|32p . When h and e 
are sufficiently small, for j = 1,2, let gj be constructed according to \39\l with the 
CGO solutions Uj . dj and dj are two sets of observations of the internal data on 
0. 

Restricting to Q, dj = dj implies that (fi,q) = 

Proof. Wc have proved that, for j = 1,2, when gj is properly chosen, (i is uniquely 
reconstructed on f2, i.e., fi = jl. Directly, dj = dj also implies u := Uj = iij. By 
unique continuation, u cannot vanish on an open set in f2 different from the empty 
set. Otherwise u vanishes everywhere and this is impossible to satisfy the boundary 
conditions. Therefore, the set F = {\u\ > 0} n Cl is open and F = Cl since the 
complement of F has to be empty. By continuity, this shows that q = q on Q. □ 

Since the coefficient q in the Schrodinger equation is unknown, the CGO solu- 
tions iij, j — 1,2, cannot be explicitly determined. Therefore, although we know 
that gj can be chosen from a open set close to Uj, a more explicit characterization 
of the open set is lacking. 

Also notice that the parameters h and e need to be small to make j3 flat enough, 
while cannot be too small, otherwise g\ will be so large that the imposed illu- 
minations become physically infeasible and 172 will be so small that the imposed 
illuminations become physically undetectable. 

4 Stability of the reconstruction 

In this section, we consider the stability of the proposed reconstruction method. 
We divide the proof into two cases. When two complex- valued data are mea- 
sured, a strict convexity is assumed on the domain of interest. In another case, 
2n complexed-valued data are measured and the stability result follows without 
any geometric condition on the domain. 

4.1 Stability Result for 2 complex observations 

To avoid the case when the vector field is approximately parallel to the boundary 
<9f2+, i.e., /3(x+(x)) • n(x+(x)) is small, for a long time, we impose the convexity 
condition as established in Hypothesis (|2.2j) . 

Recall that 9 x (t) is the flow associated to f3. Assume x (t) reaches the boundary 
at x+ and at time t + , i.e., 6 x (t+) = x + G dQ. Similar notations are use for (3. We 
have the equality 

e x {t) - e x (t) = f \p{6 x {s)) - [3(e x { s ))]ds. 

Jo 
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Using the Lipschitz continuity of /3 and Gronwall's lemma, we thus deduce the 
existence of a constant C such that 

\e x (t) - e x (t)\ < ct\\p - p\\ C o ( n) 

uniformly in t knowing that all characteristics exit f2 in finite time and provided 
that 9 x (t) and 9 x (t) are in f2. 

Such estimates are stable with respect to modifications in the initial conditions. 
Let us define W(t) = D x 9 x {t). Then W solves the equation W — D X /3(9 X ) with 
W(Q) — I and by using Gronwall's lemma once more, we deduce that 

\W(t) - W(t)\ < Ct\\D x f3 - Dj\\ c o {nh 

for all times provided that 9 x (t) and 9 x (t) are in f2. As a consequence, since fi and 
$ are in C (0), then we obtain similarly that: 

\D%(t) - D%(t)\ < Ct\\p - p\\ cKQ)l (45) 

and this again for all times 9 x (t) and 9 x (t) are in f2. To simplify the notation, we 
define 6 = \\0 - P\\ C "(n)- 

Recall that by the definition of fl in (|32p . Let x € Cl and y £ d£l+ n d£l, then 
/3(x) ■ n(jj) > r), for some constant i] > 0, where j3 satisfies (|3T[) and n(y) is the unit 
outer normal. 

Lemma 4.1. Lei fc > 1 and assume that j3 and j3 are C (0) vector fields that are 
sufficiently flat, i.e., h is sufficiently small. Let us assume that dfl+ is sufficiently 
convex so that Hypothesis \ 2.2\ holds for some R < oo. Then, restricting to Q, we 
have that 

\\ x + - x+\\ c x(n) + 11*+ - i+\\ C k(h) ^ C WP ~ ^llc*(n)' ( 46 ) 
where C is a constant depending on h and R. 

Proof. Let x E fl. Without loss of generality, we assume t + (x) < t + (x). Draw a 
circle with radius R and tangent to Q at x+. We impose a coordinate system by 
choosing the circle center to be the origin. Equation (|37|) and (|45|) gives that 

\x+ - [9 x (t+) + P(9 x (t + ))(t + - t+)]\ < C(t + - t+)S. 

Thus, 

\9 x (t+) + P(9 x (t+))(t+ - t+)\ <R + C(t+ - t+)S. 
Directly, we have 

\x + + P(9 x (t + ))(t + -t+)\ 2 -R 2 

= \9 x (t + ) + f3(9 x (t + ))(i + - t+) - (9 x (t + ) - x + )\ 2 - R 2 

= K[(R + C(t+ - t+)6 + \9 x (t+) - x+\) 2 - R 2 ]. (47) 

where < K < 1 is defined by ratio and K € C k when \9 x (t + ) — x + \ > 0. On the 
other hand, 

\x+ + P(9 x (t + ))(t + -t+)\ 2 -R 2 

= 2R~p{6 x {t+)) ■ n{x+)(i+ - t+) + \p(9 x {t+))\ 2 (t+ - t+) 2 . (48) 
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Substituting (g7|) into (g5J|, we have 

A(t + - t + f + B{t+ -t+) + C = 0, 

where 

A=0(6 x (t + ))\ 2 -KC 2 6 2 , 

B = 2RP(9 x (t+)) ■ n(x+) - 2KCS(R + \8 x (t+) - x+\), and 
C = -2KR\e x (t+) -x+\- K\9 x (t+) - x+\ 2 . 

By the quadratic formula, 



—B + VB 2 - AAC -2C , . 

t+~t+ = ^— ; = ; (49 

2A B + slB 2 - 4AC 

Notice that when x £ (l and x+ £ dfi+ n Sf2, for h sufficiently small, j3(9 x {t+)) • 
n(x+) > r\ for some fixed rj > 0. When S is small, we see that A > and B > rj > 
for some rj . Thus, (|4"5"]) and (|49p implies 

«+ - *+ < C'\9 x (t+) -x+\< C"\\(3 - PW^y (50) 

By taking derivatives of (|4"9")l , 

fe 

^(t+ - t+) =J2<*idi\Ut+) ~ (51) 

i=0 

where |aj| is bounded when <5 is small. 

We next need to estimate d x \9 x {t + ) — x+\. When i = 0, it is directly from (|4"5j) ; 
When i = 1, we have 

9 x |0 x (i+)-o; + | < |5 x (fl a[ (t + ) - a: + )| 

= \(W(t + ) - W{t + )) + (MAU)) - m(t+)))d x t + \ 
<C||/3-£|| cl(fi)) 

where the inequality follows from (|45p and the continuity of (3. Similar estimate 
works for high order derivatives. Thus, 

P+-t+|| i ( n) < WP-PWcHny (52) 
Furthermore, the estimate of \x+ — x+| follows easily. 

\d x (x+ - x+)\ = \d x (0 x (t+) - 9 x (t+))\ 

= \(W(t + ) - W(t + )) + 0(t + )dj+ - p{t+)d x t + )\ 

<C\\p-p\\c^Siy 

where the inequality follows from (|4"5|) . (|5l2"j) and the continuity of (3 and W. High 
order derivatives follow similarly. This finishes the proof. □ 



12 



Proposition 4.2. Let k > 1. £e£ /i and // be .solutions to J^g| ) corresponding to 
coefficients (/?, 7) and (/3,7), respectively, where \37\ ) holds for both (3 and f5. 

Let us define fio = and fiQ = fi\dQ,, thus /io,Ao € C (9f2). FFe afeo assume 
h is sufficiently small and f2 is sufficiently convex that Hvvothesis \2. °A holds for some 
R < 00. Then there is a constant C such that 



1 1 A* l 1 \\c k - 1 (n) 

< C||^o||c"=(ao + ) CII/3 - P\\c"{h) + h - y\\c"-Hh)) + ^IIMo - Mo||c=(9n+)- (53) 

Proof. By the method of characteristics, n(x) is determined explicitly in (|43[) . while 
/i(x) has a similar expression. 



\fj,(x) - p.(x)\ < 



+ 



e -/o + W 7(^(s))d S 



Applying Lemma |4.1[ we deduce that 

|L»^ 1 [Ai (a;+(a;))-Ao(i+(a;))]| < ||Mo-Mo||c^-i(9n)+ c 'llA i o||c'=-i(^+)l!/ 3 -/5|lc fc - 1 (n)- 

(54) 

This proves the no(x+(x)) is stable. To consider the second term, by the Leibniz 
rule it is sufficient to prove the stability result for Jq + j(6 x (s))ds. 

Assume without loss of generality that t + (x) < t + (x). Then we have, applying 

(USD, 

ft+{x) rt+(x) 



[ 7 (6 x (s))-j0 x (s))]ds = 



[(7(0*00) - 104s))) + (7 - l)0 x (s))]ds 



< c h\\c°(h)\\P - P\\c°(h) + c h - 7ll c o(n)- 

Derivatives of order k — 1 of the above expression are uniformly bounded since 
t+(x) <G C k (Cl), 7 has C k derivatives bounded on f2 and x (t) is stable as in (|45j) . 
It remains to handle the term 



v(x) 



t+(x) 



j{e x ( s ))d s . 



t+(x) 



(3 and 7 are of class C k (il), then so is the function x —> ; y(9(s)). Derivatives of 
order k — 1 of v(x) involve terms of size t+(x) — t+(x) and terms of form 



D^it+D*- 1 



^(6 x {t + ))-t + D k x -^{e x {t+))), <m< k-1. 
Since the function has k — 1 derivatives that are Lipschitz continuous, we thus have 

\D k - 1 v(x)\<C\\i+-t + \\ ck - lm . 
The rest of the proof follows Lemma |4~T] □ 
With all prepared, we are ready to prove the following stability result. 
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Theorem 4.3. Let k > 3. Assume that (fM,q) and (p,q) are in V and that \gj — 
Uj\dn, j = 1; 2, are sufficiently small. Let /i, ft be solutions of with (/3, 7), ([3, 7) 
and /i sufficiently small. Assume on the boundary ,fiQ = fio- Then we have that 

+ lk-5llc*-3(fj) ^ ^ll(CMn)) 2 - ( 55 ) 

Proof. The first part follows directly from (|35p and Proposition 14.21 This provides 
a stability result for v = and thus for itj = i^dj . Notice uj is non- vanishing on 
Q. So when choosing \gj — Uj\dn sufficiently small, the arguments in Q39p and (|40[) 
show that is non- vanishing in Q. Thus Auj + ujq = gives the stability control 
of g. □ 



4.2 Stability Result for 2n complex observations 

We now consider the case that An real-valued observations are taken, with which 
we can construct 2n sets of complex-valued boundary and internal data, denoted 
as 9i 2 = {di^A} an d d{ 2 = {di,^}. For the rest of this section, we choose 

j = l,...,n. 

Let H be a hyperplane in K n \ch(0). Choose x 3 € H, such that {x 3 — a; 1 } form 
a basis of thus spanja/' — a; 1 } has dimension n — 1. Then for Va; € fi, {a:-' — a;} 
form a basis of R™. In fact, since x ^ H, 

spaxi{x J — x} = span{a;-' — x , x 1 — x) 

has dimension m. 
Define 

fx — x 3 \ 

ip> = loglx-a^l, ^ = dist | -,w , j = l,...,n. 

\\x-x 3 \ J 

Then the matrix B v := (\7ip 3 ) is invertiblc. 

Corresponding to x 3 , tp 3 and V' 3 1 we can define the front and back sides of the 
boundary, dfl J + and dfl 3 _, by ©, and define the CGO solutions, u J 12 , by (f3"U)) . 

Let us now choose boundary conditions (7^ 2 close to 2 on dCl, precisely, by 
(j3"9")l . With internal data d{ 2 , we can define (3 3 by ([55]) and (j3"6")l . 

Proposition 13.41 shows that the matrix i? = {(3 3 ) is close to an invertible matrix 
B v on a subregion f2. Therefore, we can make B invertible on Q with inverse of 
class C k (il) by choosing h sufficiently small. Consequently, (|42l) can be rewrite 

v> + A^ = 0, (56) 

where A is a vector-valued function in (C fc (fi)) n . Finally, the construction of A is 
stable under small perturbations in the data d 3 . Indeed, let T and T be two vector 
fields constructed from the internal measurements d{ 2 and d{ 2 , respectively. Then 
when h is sufficiently small, 

ll r - r!l(c fc (o))™ < \\d- d\\( C k+i(n^n- (57) 

Now we consider equation (|56[) with boundary condition fi = fiQ. Assume fl is 
bounded and connected and <9f2 is smooth. Let x € fi. Find a smooth curve from 
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a; to a point on the boundary. Restricted to this curve, (fijo]) is a stable ordinary 
differential equation. Keep the curve fixed. Let /z, jl be solutions to (|56[) with respect 
to r, r, respectively. Assume fiQ = (1q on dfl. By solving the equation explicitly 
and applying (|5"7|) . we find that 

frWcnn) ^ ll d -^llc+i(n)- ( 58 ) 
We can now state the main theorem of this section. 

Theorem 4.4. Let k > 2. Assume that (fj,,q) and (fi,q) are in V and that we have 
In well-chosen complex valued measurement, such that \g{ 2 — u\ 2 |an> j — 1, • ■ • , n, 
are sufficiently small. Let [i, ji be solutions of \ 56}) with T, T and h sufficiently small. 
Assume on the boundary, fj,Q = Mo- Then we have that 

1 1 A* - Mllc fc (o) + Ik - q\\c k -2{n) ^ C \\ d - ^(C+H")) 2 - ( 59 ) 

Proof. The first result is directly from (|58|) . The proof of the stability of q is exactly 
that same as in the proof of Theorem 14.31 □ 



5 Proof of Main results 

In this section, we return to the inverse diffusion problem with internal data ([T]) and 
prove the main theorems. We will invert the standard Liouville change of variables 
mentioned at the beginning of section [31 Recall that by defining v = \f~Du, we had 
PfTOllTTl). which induce 

- AVD - qV15 = fi. (60) 

Section [3] allows us to reconstruct fi and q, while D is given on dQ. Thus we can 
solve for \[~D from (|60|) and then a a = fiy/T). 

The uniqueness of thus a solution \[D of (|60|) is based on that is not an 
eigenvalue of A + q. It is enough to prove that (D,a a ) G M. implies (q, /i) £ V . 
Indeed, the inverse of A + q is compact and \[D e Y = H^ +k+2+t (fl) C C k+2 (Q). 
These imply A + q is surjective and thus does not have as an eigenvalue by the 
Fredholm alternative. 

So far we proved the unique reconstruction of D,o~ a form internal data for well- 
chosen partial boundary illuminations as stated in Theorem 12.11 

Proof of Theorem \2.3[ Since a = \i\f~D, we only need to prove the stability of \f~D. 
When k > 3, we have the stability of the reconstructions of q € C k and [i e C k+1 in 
Theorem 231 while the boundary illuminations are chosen from a open set in C k,a 
with a > 1/2. According to (|60"]) . we have 

-(A + «)("/D- V r D)=(l-fl+(q-q)v r D. (61) 

By elliptic regularity, we deduce that (y/D— \/~D) is bounded in C k (ti), and hence 
the theorem. □ 

The proof of Theorem 12.41 will be the same as above. 
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